Hierarchical mesophases of vortex matter in layered superconductors 
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We demonstrate the possible formation of hierarchical mesophases of vortex matter in layered 
superconducting systems where there are variations of interlayer thicknesses and layers are made 
of different superconducting materials. We show that because inter-vortex forces feature multiple 
length scales in this case the magnetic response features the formation of mesophases such as clusters 
of vortex clusters, concentric vortex rings, vortex clusters in a ring, and vortex stripes in a cluster. 
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Condensed matter physics has long been concerned 
with explaining phenomena that result from competing 
interactions, covering a wide variety of topics from soft 
condensed matter systems to magnetism and ultracold 
atoms (for a recent overview see e.g. [1, 2]). The richest 
pattern forming systems are those with several length 
scales. For example, structure formation in systems with 
two-scale repulsive interactions is highly relevant in hard 
condensed matter systems [3, 4], nuclear matter [5], and 
in colloids and other soft condensed matter systems [6,7]. 

By contrast research on the magnetic response of type- 
2 superconductors traditionally deals with structure for- 
mation of vortices [8] . Although the interaction between 
vortices in these materials has a simple monotonically 
repulsive form, the vortex matter exhibits a plethora of 
interesting phase transitions and structure formation [9] . 
Moreover, the vortex states, especially in the presence 
of pinning, are critically important for technological ap- 
plications of superconductors, where control over vortex 
matter in many cases amounts to control of dissipation. 

Recently, the possibility of more complicated inter- 
vortex interactions in newly discovered systems has at- 
tracted much attention in various contexts: in multi- 
component superconductors [1, 10-26], superfluids [27], 
vortices in dense nuclear matter in neutron stars [28] , and 
quantum Hall systems [4] . Recently a problem which at- 
tracted interest was the phase separation in vortex mat- 
ter with long-range attractive and short-range repulsive 
inter- vortex interactions. Such forces in multi-component 
superconductors originate in the regime where there 
are several superconducting components originating from 
different bands. This gives rise to two coherence lengths, 
^1 and ^2 , and the magnetic field penetration length falls 
between them: < A < ^2 [10]. This phenomena, re- 
cently referred to as "type-1.5 superconductivity" [11], is 
the subject of a recent review [29] and has stoked consid- 
erable experimental interest in pursuing a realization of 
such regimes using artificial structures made of alternat- 
ing layers of type-1 and type-2 materials. 

Here we propose that superconducting systems can 
have vortex states with several length scales of repul- 
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FIG. 1: (Color online) (a) Sketch of a vortex showing inter- 
action length scales, an is the hard-core radius and ai (0-2) is 
the inner (outer) soft-core radius, (b) Pair potential V{r) as 
a function of the interparticle separation r for a step potential 
and an analogous smooth potential, (c) Schematic drawing of 
the field lines of a vortex in a layered superconductor. Shaded 
(white) regions are superconducting (insulation) of varying 
thickness. By controlling thickness of the of insulating layers 
one can control intervortex interaction. Thicker insulating 
layers cause a wider spread of the magnetic field lines result- 
ing in the presence of an additional repulsive length in the 
intervortex interaction scales. 



sive (also in some cases attractive) interactions, where 
more complicated interaction potentials can be realized. 
Such intervortex forces should arise in layered structures 
made of combinations of type-2 and type-1 superconduc- 
tors where the magnetic field penetration length varies 
in different layers or where there are insulating layers of 
varying thickness. Also, for quasi-two-dimensional sys- 
tems an additional repulsive interaction is present due to 
the interaction of stray fields outside the sample [30] Here 
we consider the regime where vortex line tension is large 
and temperature is small so the vortices do not bend and 
we have effective two-dimensional system where vortices 
can be described by their positions in xy-plane. 

With this physical realization in mind, let us con- 
sider the simplest potential with several length scales, a 
hard sphere model with multiple shoulders. Note that 
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FIG. 2: (Color online) Final particle configuration for 
(a) three-step potential of Fig. 1(b) and (b) smoothed po- 
tential with A'^ = 2000 and p = 0.51. The corresponding 
radial distribution functions are plotted versus particle sepa- 
ration (in units of the characteristic simulation length A [34]) 
in panels (c) and (d), respectively. 



"core-softened" potentials with a single shoulder have 
been studied intensively, revealing a myriad of density- 
modulated ground states [6, 7, 31-33]. In Figs. 1(a) and 
1(b) we illustrate a particle with an impenetrable hard- 
core radius ah and two repulsive shoulders at r = (Ti and 
CT2 with heights ei and £2, respectively. 

We show the final configuration from Monte Carlo 
(MC) simulations [34, 35] of our hard-sphere model in 
Fig. 2(a) for Ny = 2000 particles and density p = 0.51. 
Here the system forms a mesophase: namely the particles 
order on three different length scales: (1) the particles 
form a tightly bound cluster, (2) the clusters are them- 
selves bound into a conglomerate structure (hereafter re- 
ferred to as a superclustcr) , and (3) the structures form 
a lattice. To analyze the underlying structure of this 
phase, we show the radial distribution function (RDF) 
g{r) [34, 36] in Fig. 2(c). The first feature in g{r) is a 
very strong peak corresponding to the nearest-neighbor 
distance inside each cluster. Because the clusters have 
such a small radius, g{r) shows only a small peak at ap- 
proximately double the nearest-neighbor distance. The 
next pronounced peak is the inter-cluster distance inside 
a superclustcr, with the subsequent peaks describing the 
distance between clusters in different superclusters. 

A more physically accurate potential is achieved by 
smoothing out the steps in the previous potential [see 
Fig. 1(b)]. In Fig. 2(b), we show the final MC config- 
uration for the smooth potential, which is a superclus- 
ter of higher radial symmetry. The RDF is plotted in 
Fig. 2(d) and has nearly identical features, although the 
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FIG. 3: (Color online) Inter-vortex pair potentials used in 
this study. Panels (a), (b), and (c) each illustrate a different 
pair potential, with the top and bottom of each panel repre- 
senting a different set of parameters [34] . The insets are close 
up views of the potentials for small interparticle separation. 
The legend indicates the corresponding panels in Fig. 4. 



higher symmetry of the ground state results in smoother 
peaks. The conclusions that follow are (i) multiple re- 
pulsive length scales result in a formation of hierarchical 
structures and (ii) the precise form of the potential is of 
a lesser importance in this example: the crucial aspect is 
the existence of several length scales in the interaction. 

To classify possible mesophases we need to consider 
what sorts of vortex structures exist in systems where 
there exist competing interactions that arc dominant at 
different length scales (see Fig. 3). In the context of 
Fig. 1, such systems can be realized by adding layers of 
type-1 material or alternating layers of clean and dirty 
material, while repulsive scales are tunable by controlling 
e.g. layer thicknesses. 

In Fig. 3(a), we show two potentials that both feature 
a strong repulsive core surrounded by an area of attrac- 
tion. Outside the attractive shell, there is a repulsive 
region and an attractive long-range coupling. At very 
low vortex densities, the final MC configurations for both 
of these potentials are given in Fig. 4(a) and Fig. 4(b), 
respectively. In the first case, the particles form a clus- 
ter due to their attractive interaction, the repulsive scale 
however gives this cluster a ring shape. In the second, 
the combination of attractive and repulsive scales induces 
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FIG. 4: (Color online) Snapshot of the final vortex configurations corresponding to the potentials in (a) top panel of Fig. 3(a) 
with Nv = 1000 and density p = 0.044, (b) bottom panel of Fig. 3(a) with A''„ = 1000 and density p = 0.025, (c) top panel of 
Fig. 3(b) with = 2000 and p = 0.20, (d) bottom panel of Fig. 3(b) with = 1000 and p = 0.40, top panel of Fig. 3(c) with 
A''^ = 3000 and densities (e) p = 0.25 and (f) p — 1.00, and the bottom panel of Fig. 3(c) with A'^„ = 3000 and densities (g) 
p = 0.50, (h) p = 1.25, (i) p = 1.50, and (j) p = 2.50. For panels (a), (b), (c), and (d), the long-range attraction causes all of the 
particles to form a single object and we only show a close up view. The unlabeled panels are close-up views of panels (e) and 
(f), focusing on a single supercluster and ring, respectively. The final vortex configurations in panels (a)-(d) and panels (e-j) 
are from Monte Carlo [35] and Molecular Dynamics [37] simulations, respectively, with simulation details discussed in Ref. 34. 



clustering inside the ring (a clustered ring). 

Next, we consider two potentials where the long-range 
attraction is extremely weak and we vary the potential at 
intermediate length scales [see Fig. 3(b)]. The final MC 
configurations for these potentials at densities p ~ 0.2 
and p = 0.8 are shown in Figs. 4(c) and 4(d). In the first 
case, the particles form a single large supercluster with 
one critical difference from the case shown on Fig. 2: the 
size of the constituent small clusters is modulated by the 
distance to the center of the cluster, going from a maxi- 
mum of 4 vortices per cluster in the center to a shell of 
single vortices at the edge. In the second case, when the 
interaction in the intermediate region is modified, the lo- 
cal phase of the cluster varies with distance from the cen- 



ter: namely as one goes from the center of the cluster to 
the edge one encounters regions corresponding to vortex 
lattice, vortex stripes, and vortex voids phases. Here the 
long range attractive interaction makes the vortex den- 
sity gradually increase towards the center of the cluster 
leading to a sequence of mcsophases which optimize the 
interaction associated with repulsive short-range scales. 

The third pair of potentials [see Fig. 3(c)] we examine 
feature a moderate repulsive core surrounded by an at- 
tractive well and have a long-range repulsive interaction. 
When the well is strong, the final vortex configurations 
from Molecular Dynamics (MD) simulations [34, 37, 38] 
are circular superclusters at a density p = 0.25 or concen- 
tric rings at a density of p = 1 [illustrated in Figs. 4(e) 
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FIG. 5: (Color online) Radial distribution function g{r) cor- 
responding to the phases illustrated in Fig. 4. 

and 4(f)]. Again the short-range structure is determined 
by repulsive length scales. However, when the attrac- 
tive well is weakened significantly more regular repulsion- 
dominated vortex phases appear at higher densities [see 
Figs. 4(g)-(j)]: a triangular lattice, a pair vortex lattice, 
stripes, and voids, which is consistent with a dominance 
of short-range two-scale repulsive interactions. 

To better understand the structure of these phases, let 
us examine the RDF for each phase, which are shown 
in Fig. 5 (note that the ordering of the panels matches 
the ordering of Fig. 4). For the ring phase of Fig. 4(a), 
g{r) has three pronounced peaks indicating the nearest- 
neighbor distance, next-nearest-neighbor distance, etc. 
inside the ring. In between the peaks, g(r) remains finite 
because the vortices do not form an even lattice inside 
the ring. For the clustered ring phase [see Fig. 4(b)], the 
first peak is very pronounced and narrow, indicative of 
the particles in each cluster being roughly equidistant. 
The second peak characterizes nearest-neighbor distance 
between clusters, and the subsequent peaks the distance 
between next-nearest-neighboring clusters, etc. The long 
distances between the peaks indicate that the clusters arc 



small compared to their separation. 

The modulated supercluster of Fig. 4(c) likewise has 
a sharp, narrow peak in g(r) representing the nearly 
equidistant particle separation inside each cluster. For 
the stripe/ void-rich phase [see Fig. 4(d)], we observe sev- 
eral prominent peaks in g(r) that are largely consistent 
with the RDF for a triangular lattice, with broadening 
of the peaks due to the mixing of phases. 

The superclusters shown in Fig. 4(e) possess short 
range periodicity inside each cluster, resulting in two nar- 
row peaks in g{r) which are so close together that the 
second peak appears as a shoulder. The successive peaks 
illustrate the distance between clusters and are broad- 
ened because of the finite size of each cluster. Because 
the clusters are widely separated, there are additional 
broad peaks in g{r) that occur for large r and describe 
the supercluster separation. The RDF of the concentric 
ring phase [pictured in Fig. 4(f)] is remarkably similar 
to the g{r) for the supercluster phase [Fig. 4(e)]. Unlike 
the supercluster phase, g{r) remains finite for r smaller 
than the diameter of the outer ring due to the particles 
spreading out evenly throughout each ring. 

Finally, we discuss the radial distribution functions of 
more conventionally ordered phases. The lattice phase 
[see Fig. 5(g)] possesses much more long-range order than 
all other phases considered, with peaks at r = a, ^/3a, 
2a, \/7a, 3a, . . ., where a is the nearest-neighbor dis- 
tance for a triangular lattice. The dimer lattice phase [see 
Fig. 5(h)] has a clearly defined peak describing the dimer 
separation. Because the pairs prefer to line up end-to- 
end, the peaks that would describe the triangular lattice 
are broadened significantly and the long-range order can- 
not be observed in g{r). The stripe phase [see Fig. 5(i)] 
has several regularly spaced peaks coinciding perfectly 
with the separation of particles along each stripe. Here, 
broadening occurs due to both bending of the stripes 
and the presence of other stripes. The void phase [see 
Fig. 5(j)] only has short-range periodicity, as evidenced 
by the pronounced features for small r. 

In this paper, we studied vortex states in general 
layered superconductor-insulator-superconducting struc- 
tures, made of different superconducting layers. The 
vortices are subject to interactions with multiple length 
scales. We have shown that these layered systems have 
an unusual magnetic response: vortex supercluster struc- 
tures, which can consist of clusters of clusters, rings, clus- 
ters in a ring, or have coexistence of stripes, voids, and 
lattice phases. This can provide an experimental tool to 
deduce information about vortex interactions from ob- 
servation of their ordering in real experiments. Besides 
that it indicates that one can use layered superconduct- 
ing structures for the realization and study of rich and 
unique pattern forming systems. In further studies we 
plan to address effect of disorder and pinning on these 
vortex states, as well as their manifestations in the trans- 
port properties of layered structures. 
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